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FINITENESS THEOREMS AND ALGORITHMS FOR 
PERMUTATION INVARIANT CHAINS OF LAURENT LATTICE IDEALS 

CHRISTOPHER J. HILLAR AND ABRAHAM MARTIN DEL CAMPO 

Abstract. We study chains of lattice ideals that are invariant under a symmetric group action. In 
our setting, the ambient rings for these ideals are polynomial rings which are increasing in (Krull) 
dimension. Thus, these chains will fail to stabilize in the traditional commutative algebra sense. 
However, we prove a theorem which says that "up to the action of the group", these chains lo- 
cally stabilize. We also give an algorithm, which we have implemented in software, for explicitly 
constructing these stabilization generators for a family of Laurent toric ideals involved in applica- 
tions to algebraic statistics. We close with several open problems and conjectures arising from our 
theoretical and computational investigations. 



1. Introduction 

In commutative algebra, finiteness plays a significant role both theoretically and computationally. 
An important example is Hilbert's basis theorem, which states that any ideal / C i? in a polynomial 
ring R = C[xi, . . . ,Xn] over the complex numbers C (or more generally, over any field K) has a 
finite set of generators G = {gi, . . . ,gm}- 

I = {G)r ■.= giR + --- + g,nR- 

In other words, C[xi,... Noetherian ring. Equivalently, any ascending chain of ideals 

-^1 ^ -^2 ^ • • • in C[xi, . . . ,Xn] stabilizes (i.e., there exists an A'' such that In = In+i = •••)• 
This result has many applications in the algebraic theory of polynomial rings (e.g. the existence of 
finite resolutions |15l p. 340]), but it is also a fundamental fact underlying computational algebraic 
geometry (e.g. termination of Buchberger's algorithm in the theory of Grobner bases [10^ p. 90]). 

In many contexts, however, finiteness is observed even though Hilbert's basis theorem does not 
directly apply. A motivating example is the (non-Noetherian) ring R = C[xi,X2, . . .] of polynomials 
in an infinite number of indeterminates X = {xi, X2, ■ ■ ■}, equipped with a permutation action 
on indices. More precisely, the symmetric group Sp of all permutations of the positive integers 
P := {1,2,...} acts naturally on R via: 

(1-1) O-f {Xsi, Xsi) := f{x„(^si),...,X„(^Si)), CT£&p,feR. 
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Although many ideals in the ring R are not finitely generated, an important subclass still admit 
finite presentations. Call an ideal / permutation-invariant if it is fixed under the action of (3p: 

6p/:= M:aG6p, /£/} = /. 

It is known that for every such permutation-invariant I Q R, there is a finite set of generators 
G = {gi, . . . , gm} C / giving it a presentation of the form: 

As a simple example, the ideal M C C[xi,X2,...] of polynomials without constant term has the 
finite presentation M = {&fXi)c[xi,x2,...] even though it is not finitely generated. 

The above finiteness property for the ring C[xi,X2,...] was first discovered by Cohen in the 
context of group theory [8] (see also [9] for algorithmic aspects), but seems to have gone unnoticed 
in the commutative algebra community until its independent rediscovery recently in [2]. Gener- 
alizations and extensions of this result have since been applied to unify several finiteness results 
in algebraic statistics [SS] as well as help prove open conjectures in that field (notably, the inde- 
pendent set conjecture [261 125]. finiteneness for the A;-factor model [T2], and, more recently, that 
bounded-rank tensors are defined in bounded degree ^13j). 

In this paper, we derive new finiteness properties for certain classes of polynomial ideals that 
are invariant under a symmetric group action. Motivated by an algebraic question of Dress and 
Sturmfels in chemistry [U Section 5] , we prove that invariant chains of lattice ideals stabilize up to 
monomial localization (see Theorem 11.31 below). This general result gives evidence for Conjecture 
5.10 in [4] (stated as Conjecture 14.21 below). Moreover, for the specific chains studied there (in [U 
Section 5.1]), we present an algorithm for explicitly constructing these generators (see Theorem 11.71 
and Algorithm [T] below) . Our results also have potential implications for algebraic statistics. To 
prepare for the precise statements, however, we need to introduce some notation. 

Given a set S, let 65 denote the group of permutations of S. We shall focus our attention 
primarily on the sets S = [n] := {1,2,..., n} and 5" = P := {1, 2, . . .}, the set of positive integers. 
In these cases, we write 6„ and 6p, respectively, for the symmetric groupsQ Given a positive 
integer A; > 1, let [S]^ be the set of all ordered A;-tuples u = {ui, . . . , n^), and let (5)'^ be the subset 
of those with pairwise distinct ui, . . . ,Uk- When S = [n], we write [n]^ and (n)^ for [S]^ and {S)^, 
respectively. 

The symmetric group 65 acts on [S]'^ naturally via 
(1.2) a{ui,...,Uk) := {cr{ui),...,a{uk)), for <t G 65; 

and this action restricts to an action on (5)'^. 

Write Xs ■= {xs : s G 5} for the set of indeterminates indexed by a set 5, and let ]K[Xs] denote 
the polynomial ring with coefficients in a field K (e.g., C or M) and indeterminates Xs- The action 
of any group © on 5 induces an action on Xs, which we extend to an action on ^[X^] as in (II. ip . 

We are interested here in the highly structured G -invariant ideals of IC[X5] (simply called in- 
variant ideals below if the group (3 is understood); these are ideals / C for which GI = /0 
Guised in various forms, invariant ideals of polynomial rings arise naturally in many contexts. For 
instance, they appear in applications of polynomial algebra to chemistry [33^ HI [12], finiteness of 

-'^We embed &„ into Sm for n < m in the natural way. 

^In the language of [?], invariant ideals are also the K[Xs] * (3-submodules of K[Xs], where ]K[Xs] * 6 is the skew 
group ring associated to ]K[Xs] and 6. 
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Statistical models in algebraic statistics and toric algebra [MllMlEaillESlilEaElElIISllSnillSl 
\W\ [131 EH 123 , and the algebra of tensor rank [TBJ . 

Given an ideal / C of a polynomial ring R = M.[Xs], let denote the localization / of I 

with respect to the multiplicative set of monomials of R (including the monomial 1). In particular, 
R^ is the ring of Laurent polynomials in the indeterminates of R, and any ideal ICR lifts to an 
ideal C i?^, which we call a Laurent ideal. In simple terms, the ideal consists of elements 
of the form gh~^ where g G I and /i is a monomial of R (see e.g. [15]). An action of a group 6 of 
automorphisms that permute the indeterminates of R extends naturally to an action on R^: for 
cr G S and gh^^ € R^, we can define a{gh~^) := a{g)a{h)^^ € R^. In this way, any ©-invariant 
ideal / lifts to an S-invariant ideal C R^. As above, for a subset G Q R, we let {G)ji denote 
the ideal generated by G over R. 

In this paper, we work with localized (Laurent) ideals because they allow us to prove very general 
finiteness theorems in cases where no other known techniques are able to produce such results. 

Fix a positive integer k. In what follows, we are primarily concerned with the polynomial rings 
(and their localizations): 

(1.3) TZn := IC[X[„]fe], 7^p := K[X[p]fe] = |J 7^„; i?„ := K[X(„^.], Rp = \J Rn, 

ngP neP 

and Tn := ]K[ti, . . . , Since the set [n]^ sits naturally inside [m]'' for n < m, we have an 
embedding of rings TZn ^ Tim', similarly, Rn C R^. Our main objects of interest will be ascending 
chains Iq of ideals In ^ Hn (simply called chains below): 

(1.4) /o :=/i C ... . 

In general, a chain of ideals (jl.4p will not stabilize in the sense of Hilbert's basis theorem because 
the number of indeterminates in TZn increases with n. However, if the ideals comprising a chain are 
©-invariant, we may still be able to find an N such that all the ideals Im-, In+i, ■ ■ ■ are the same. 
We now make these notions precise (with corresponding definitions for Laurent ideals and the rings 

Rn)- 

Definition 1.1. A chain lo '■= Ii Q I2 ^ ■ ■ ■ of ideals /„ ^ TZn is an invariant chain if 

&mln Q Im, for all m > n. 
Definition 1.2. An invariant chain Iq stabilizes if there is an integer N such that 

{&mlN)nm = Im, for all m> N. 
Such an is a stabilization bound for the chain, and generators for I^ are called generators for Iq. 

In words, an invariant chain stabilizes when its fundamental structure is contained in a finite 
number of ideals comprising the chain. When k = 1, every invariant chain of ideals in -[T^^j-j^gijp 
stabilizes [H[25]. However, the corresponding fact fails to hold for k >2 (e.g., see [H Proposition 
5.2] or |25t Example 3.8]), and more refined methods are required to detect chain stabilization. 

In many applications, the invariant chains consist of toric ideals, so we shall focus our attention 
here on the slightly more general class of lattice ideals (see Section [3] for definitions). For instance, 
the independent set conjecture in algebraic statistics |26^ Conj. 4.6] concerns stabilization for a 
large family of toric chains. 
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Our first main result asserts that invariant chains of lattice ideals stabilize locally, and it is 
similar to a chain stabilization result used in a recent proof |25] of the independent set conjecture. 
We prove this result in Section [3] using ideas from order theory as described in Section [2j 

Theorem 1.3. Every invariant chain ■= Q - - - of Laurent lattice ideals C TZ^ 

(resp. C R^) stabilizes. 

Although this result is quite general, our proof is nonconstructive. In applications, however, one 
usually desires bounds on chain stabilization. Our second main result restricts to the rings Rn and 
provides a stabilization bound for the special case of Laurent toric chains induced by a monomial 
[H Section 5.2], which we study in Section HI These toric ideals appear in applications to algebraic 
statistics [171 ES] and voting theory [TTj . 

Theorem 1.4. Let f G . . . ,7/^] be a monomial of degree d in k variables. For each n > k, 

consider the (toric) map: 

4'n '■ Rn ^ ?n, ^(ni ^ f (fui ) • • • 5 tu^. ) • 

Let In = ker^n, and let 1^ be the corresponding Laurent ideal. Then N = 2d is a stabilization 
bound for the invariant chain 1^=1^'^ -^fc+i ^ ' ' ' of Laurent ideals. 

Example 1.5. Let k = 2 and suppose that / = yfy2 £ ^[yi,y2]- For every n > 2, the map is 
defined by (pn{x(ij)) = tftj for G (ra)^. Theorem 11.41 asserts that if = 2 • deg(/) = 6, then 
the generators of form a generating set for the whole chain up to the action of the symmetric 
group ©m,; that is, for all m > 6, we have {&ml^) R,n = ^m- For instance, when m > 9, we observe 
that X(3_9)X(7^9) - X(3j)j;(9^7) G Im (thus, in I^) since 

Thus, by Theorem 11.41 there exist permutations ai,...,aT G ©m, elements gi,...,gr G I^, and 
polynomials hi, . . . ,hr G R^, such that ^(3 9)^(7 g) — 3^(3 7)^(9 7) = hiUigi + • • • + hj.argr- Theorem 
[LT] below, provides a method for finding such polynomial combinations in general; in this case, one 
possibility is r = 1, /ii = 1, ai = (1 3 9 2 7) G ©m, and gi = X(^i^^)X(2,z) - x {i^2)X {z,2) ^ I& ■ For more 
details on this example (including an explicit set of generators for /^), see Section 21 □ 

Remark 1.6. Rather surprisingly, it is still an open question whether the (non-Laurent) toric 
chain /□ stabilizes in Example 11.51 and more generally, for any monomial / that is not square- free. 
Section [6] discusses more open problems of this nature. 

In the development of the proof of Theorem 11.41 we also found an algorithm for computing these 
generators. 

Theorem 1.7 (Algorithm (ll) . There is an effective algorithm to compute a finite set of generators 
for the Laurent chains in Theorem \1.4\ 

The first step of the algorithm in Theorem 11.71 is to embed a toric ideal into a Veronese ideal 
in a larger polynomial ring and use the fact that the latter is generated by quadratic binomials. 
A second procedure replaces the extra indeterminates of the larger ring by special quotients of 
monomials involving only indeterminates of the original polynomial ring. In turn, this reduces to 
an integer programming problem, which we solve explicitly. The following example illustrates some 
of the main ideas involved. 
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Example 1.8. (Continuing Example II .Sp . Consider the polynomial rings R'^ := -Rn [3^(1,2,3)] ™ 
extra indeterminate 2,3)) ^iid extend 0„ to a map 4>'^ : R'^ — )• T„ by setting i;^^(2;(i,2,3)) = ^1^2^- 
Notice that \i h ^ 1^ then h € ker(/>J^, and also that 

'?^n (2^(1,2,3)) = </'n(a^{l,3)a:^{2,3)) = <P'n{x{l,2)X{za)) = ^¥¥1- 

Thus, pi := 2;(i^3)2;(2,3) — x^-^ 2 3) ^^i"^ P2 := 3;(i^2)2^{3,2) ~ ^^(i 2 3) in ker 0^ (for n > 3). Consider 
any generating set for ker0^ which contains pi,P2', then, each g ^ In can be expressed in terms of 
these generators. For instance, 

a = X{l^z)X(2,2,) - a;(l,2)2;(3,2) = (2^(1,3)2^(2,3) " 2^^.2,3)) " (2^(1,2)2;(3,2) " 2^^1,2,3)) ^ ker(/>^. 

Next, notice that 

n r;N ^' ( \ + + + Mx(l,2))Mx{'i,l)) , / 2:(i,2)2^(3,l) 

(1-5 ^^n (2^(1,2,3) = ^1*2*3 = —, ^ = 

<Pn(2^(l,3)) V 2;(i,3) 

Therefore, if we replace X(i^2,3) by ^'-^'^'i^^^'^^ in the two generators pi and p2 above, we obtain 

two elements pi, P2 € In which also generate g. More generally, if we can find a finite set of 
generators for the chain of ideals ker 4>'^, then we would have generators for the chain of ideals I„ 
up to monomial inversion. 

Identity (jl.Sp was discovered by solving the following integer programming problem (described 
more fully in Example 14. Sp . The exponent vector of tit2t^ \s u = (1, 1, 1, 0, . . . , 0) S and for any 
(i,j) G {n)'^, the exponent vector of 4>n{x(^i^j^) = tftj is 

Wi^j := (0,. .. ,0,2,0,. ..,0, 1,0,. .. ,0) gZ", 

in which the nonzero components of Wij are the ith and j'th with respective values 2 and 1. To find 
an expression such as (jl.Sp . we needed to write u as an integer linear combination of the vectors 
Wi^j (this is done in general in Lemma l4.4p . □ 

The most recent finiteness result along the lines of Theorems 11.31 and 11.41 can be found in the 
work of Draisma and Kuttler [13]. There, they prove set-theoretically that for any fixed positive 
integer r, there exists d G N such that for all p G N, the set of p-tensors (elements of Fi • • • (8) Vp, 
where each is a finite dimensional K-vector space) of border rank at most r are defined by the 
vanishing of finitely many polynomials of degree at most d (when r = 1 these polynomials define 
toric ideals). The authors of [T3] also realized the ideals defined by these polynomial equations as 
invariant chains under the action of the semi-direct product of 6p with the general linear group 
GL{VY , and they conjectured [131 Conjecture 7.3] stabilization. The case r = 1 was proved by 
Snowden in [36]. The results of |13j extend those of Landsberg and Manivel in [29], where they 
show set-theoretically that p-tensors of rank at most 2 are defined by polynomials of degree 3 (the 
(3 X 3)-subdeterminants of all the flattenings) regardless of the dimension of the tensor. We note 
that an ideal-theoretic proof of this last fact was recently discovered by Raicu [32] . 

While the general problem of deciding which chains of ideals stabilize seems difficult, it is possible 
that every invariant chain of (non-Laurent) lattice or toric ideals stabilizes, and Theorem 11.31 
provides evidence. However, even for the special case studied here of a toric chain induced by 
a monomial, this is not known [U Conjecture 5.10] and appears to be a difficult problem (although 
it is true for square-free monomials [T, Theorem 5.7]). We pose the following open question. 
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Problem 1.9. Does every invariant chain of lattice ideals (resp. toric ideals) stabilize? 

The outline of this paper is as fohows. In Section [21 we introduce the order theory required for 
proving Laurent lattice stabilization (Theorem II. 3p in Section [3l Next, Section [5] contains a proof 
of Theorem 11.41 using some ideas from toric algebra and integer programming. Also found there 
is another approach to constructing Laurent chain generators in Theorem 11.41 (e.g.. the generators 
alluded to in Example ll.Sp which can produce smaller generating sets than those found by Algorithm 
[TJ Section [5] contains a discussion of Theorem 11.71 and Algorithm [TJ Finally, in Section [6] we present 
some open problems and conjectures arising from our computational investigations. 

2. Nice Orderings 

In this section, we explain the ideas from the theory of partial orderings that are needed to prove 
Theorem 11.31 A well-partial- ordering < on a set S is a partial order such that (i) there are no 
infinite antichains and (ii) there are no infinite strictly decreasing sequences. One can check that 
this naturally generalizes the notion of "well-ordering" to orders < which are not total. 

Let S be a group acting on a set S (a &-set), and suppose that < is a well-ordering of S. For 
s G S and o" G (3, let s< := {t G S" : t < s} and (Ts< := {at : t < s}. We define a partial ordering ^ 
on S as follows: 

(2.1) s ^ t :<^=^ s <t and there exists o" G (3 such that as = t and (TS< C t^. 

A group element a G & verifying (12. ip is called a witness of the relation s ^ t. An example of this 
construction can be found in Example 12.11 

Call the well-ordering < of S" a nice ordering if ^ is a well-partial-ordering. Many naturally 
occurring S-sets have nice orderings. For instance, the set of fc-element subsets of P with the 
natural action of (3 = (3p has a nice ordering [2]. Camina and Evans studied the ring-theoretic 
consequences of nice orderings in [7J, inspired by the ideas in [2]. They showed that if 5 has a 
nice ordering, then the ]K[(3] -module M.S is Noetherian over the group ring ]K[S] for any field IC [71 
Theorem 2.4]. We shall prove that \P]^ also has a nice ordering; however, our application (Theorem 
11. Sp requires a more refined version of this statement. This refinement is given by Theorem 12.101 
below. Before proving this theorem, we first define a nice ordering of [F]^ with special properties. 

Consider Sp acting on [Ff as described in (11. 2p . We first give a total well-ordering <diex on 
[F]^ as follows. Given w = {wi, . . . ,Wk) G [P]'^, set \w\oo := msLK{wi, . . . ,Wk}- Define the degree 
lexicographic total ordering on [F]^ by 

(2.2) V <di ex 'W \v\oc ^ |wj[oo or ["ujoo — ["U^Ioo 

and V <iex w. 

Here, </ex is the natural lexicographic ordering of elements of [F]^ given by {ui,...,Uh) <iex 
{wi, . . . , Wk) :<^=^ ui = wi, . . . , Ur-i = Wr-1 and Ur < Wr for some r G [fe]. 

Notice that for every w G [P]'^ there are only finitely many v G [P]'^ such that v <diex w; hence, 
<diex is a well-ordering of [P]'^. The well-ordering <diex induces the partial order :<diex as in (j2.ip . 

Example 2.1. With the above definition of <diex for [P]^, we have (2, 3) <diex (2, 4) and (2, 3) <diex 
(3,1). Moreover, when [P]^ is equipped with the action of Sp, we claim that (2,3) :<diex (2,4). 
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Represent the elements of Sp in cyclic notation so that (3 4) • (2,3) = (2,4), and observe that 
(3 4).(2,3)<,,^^ = (34) .{(1,3), (2,1), (2, 2), (1,2)} 
= {(1,4),(2,1),(2,2),(1,2)} 

C {(1, 1), (1, 4), (3, 1), (1, 3), (2, 3), (3, 2), (3, 3), (2, 1), (2, 2), (1, 2)} 
= (2,4)<,,,. 

On the other hand, we have (2, 3) -:^diex (3, 1). To see this, let a € Sp be such that a ■ (2, 3) = (3, 1); 
thus a{\) > 2. Notice that for (2,1) G (2,3)<^,^^, we have a ■ (2,1) = (cj(2), a(l)) = (3,cj(1)). It 
follows that (3,2) <diex c ■ (2, 1)- Since 

(3, 1)<,,.. = {(2, 3), (1, 3), (2, 1), (1, 2), (1, 1), (2, 2)}, 

we see that (3, 2) ^ (3, 1)<,,^^ ; therefore, a ■ (2, 3)<,;^^ ^ (3, 1)<,,^^ . □ 

Although not needed for our main result, a solution to the following problem would likely be 
useful in converting the methods of this section into computational tools. 

Problem 2.2. Give a computationally efficient criteria for determining ifu :^diex v foru,v € {F]^ . 

One may also ask the following open-ended problem. 
Problem 2.3. Let S be an G-set. Characterize those total well-orderings < which are nice. 

We are now in position to show that the ordering <diex is nice. 

Proposition 2.4. The ordering :<diex of is a well-partial-ordering. 

The proof of this proposition uses a special case of a result of Higman [23\ 131). which we state in 
the following lemma. Recall that a strictly increasing map (p : [m] — )■ [n] satisfies ip{i) < (p{i + 1) 
for all i. 

Lemma 2.5 (|23j). Let S be a finite set. The following ordering <h on the set S* of all finite 
sequences of elements of S is a well-partial-ordering: 

{3 tp: [m] — [n] such that ip is 
strictly increasing and Xi = y^(^i) 
for all i G [m] 

Proof of Proposition \2.4\ Let S := {0, 1, ... , k}. First order [F]^ x S* by the product of the or- 
derings <diex and <h on [F]^ and S*, respectively. Then [F]^ x S* is well-partial-ordered, by 
Higman's Lemma (the product ordering of two well-partial-orderings is a well-partial ordering). 
For w = {wi, . . . ,Wk) G [IP]'^, set n := \w\oo', also, let w* := {wl, . . . ,w^) G S* be given by 

w* ■= j for i = 1, . . . ,n. 

Wj=i 

To prove that :<diex is a well-partial ordering on [F]^, it suffices to show that the map w i— >■ 
{w, w*) : [F]^ — > [F]^ X S* is an order- embedding; that is, if v <diex w and v* <h w*, then v :<diex 
for all v,w e [F]''. 

Suppose that v <diex w and v* <h w*, and let 

^ — l^looi ^ — l^loo? tliGn tliGrc exists a- function 
(p: [m] — > [n] strictly increasing such that = for all 1 < z < m. Since (p is injective, it can 
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be extended to a permutation a € Sp. We claim that v :<diex w via witness a so that av = w and 



We first verify that av = w. For i £ {1, . . . , k}, let I = Vi < m. Notice that w*^^-^ = v*, and so 
together with the definition of v* , we have 

v*i =i + ^j ^ K[i) =i+ j- 

Vj=l Wj=a{l) 

In particular, Wi = thus, a{vi) = a{l) = Wi and so av = w. 

Now, suppose u <diex V. Since a and ip agree on {vi, . . . -.Vk}-, it follows that |(Tti|oo < Ic^loo = 
|w|oo = n. To show au <diex w, it suffices to verify this when |ti|oo = \v\ao-, as the other case 
follows from ip being strictly increasing. If |n[oo = m = \v\oo and u <diex v, there is an r G [k] 
such that ui = vi, . . . ,Ur-i = Vr-i and Ur < Vr- Therefore, a{ui) = wi, . . . ,a{ur^i) = Wj—i 
and a{ur) < a{vr) = Wr as a is strictly increasing. Thus, au <diex w and so av<^^^^^ C w<^^^^^ as 
required. □ 

Remark 2.6. Higman's lemma is also a key element in all known proofs of the finiteness result for 
6p-invariant ideals of C[xi, X2, . . .] that was mentioned in the introduction. 

The following result also follows from the proof of Proposition 12.41 

Corollary 2.7. The ordering :<dlex of (F)^ is a well-partial- ordering. 

Proof. The same proof as Proposition 12.41 works just by noticing that, in this case, w* = j if Wj = i 
or otherwise. □ 

Not all natural orders are nice as the following example demonstrates. 

Example 2.8. Define the reverse lexicographic ordering <reviex 

on as follows: 

(2.3) {Ui,. . . ,Uk) <revlex (Wi,- ■ ■ ,Wk) Uk = Wk, ■ ■ ■ ,Uk-r = W^-r 

and Wk-r~i < Uk-r-1 , 

for some r G [k]. In contrast to Proposition 12.41 the partial order <reviex is not nice. For instance, 
we have in [P]-^ the following infinite strictly decreasing sequence: 

^revlex (^i 3) '^revlex (5) 3) '^revlex i.^:'^^- ^ 

The nice ordering :<diex is useful theoretically because of the following property. 



Lemma 2.9. Let :<diex be the well-partial-ordering (12. ip induced by the nice ordering <diex of 
[P]*=. Also, let s, t G [Ff satisfy s ^diex t and |t|oo < M for some M G {0, 1, . . .}. Then there is a 
a G &M witnessing s <diex t. 



oo • 



Proof. Since s <diex t, there exists r G Sp such that ts = t and ts<_^i^^ C t<^,^^. Let M = \t 
Construct a G ©m by setting a{i) := T{i) if t(z) < M and then extending a to a permutation 
of [M]. We claim that as = t and as<:^^^^ C i<^,g^. Since s <diex t, we have |s|oo < \t\oo = M. 
Therefore writing s = (si, . . . , Sk) G [F]^, it follows that T{si) < M for each i; thus, a{s) = t{s) = t. 
Notice also that t{w) < M for all vj G s<^;g^ because for all vj G s<rf;g^, we have w <diex s which 
implies \vj\oo < |s|oo < M, and the same holds for all u G t^^i^^^. Therefore, r(wi) < \t\oo = M for 
all w G s<d;e^ and each i = 1, . . . ,k. Thus, a{w) = t{w); therefore, as<^^^^^ C t<^,^^. □ 
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If A is a commutative ring and S an S-set, we let ^45 denote the free ^-module with basis S. 
Also, let A[&] be the (left) group ring (whose elements are formal linear combinations of elements 
in & with coefficients in A [?]). The natural linear action of A[<3] on AS makes it into an 
module. The following is the refinement of the Noetherianity result from [7J that we will use to 
prove Theorem 11.31 

Theorem 2.10. Let A be a Noetherian commutative ring. For every A[&p]-submodule B C 74[P]^, 
there exists a finite set G B such that 

t 

f^B^A\vn}^ <^ ^ax,. ■ ■ ,(Ji^&rn\ gx-,- ■ ■ -.gi^G; ax,. . . ,ae,^ A with f = y^g^g-^gj. 

1=1 

Proof. Let :^diex be the well-partial-ordering of \F]^ (by Proposition 12. 4p induced by the total well- 
order ^diex from (|2.2p . A final segment of the partial order z^diex is a set F C[Ff such that ue F 
and u :<diex v implies that v F. A well-known characterization of well-partial-orderings (see e.g. 
|27j ) is that final segments are finitely generated. That is, for every final segment F, there is a 
finite set T Q F such that F = {v : 3u € T with u :<diex v}. 

If / € j4[P]'^, we define the head of /, Head(/), to be the largest nonzero element in \P]^ (with 
respect to <diex) in the support of / (those elements of [P]'^ occurring in / with nonzero coefficient). 

For the A[6p]-submodule let J C A be the ideal generated by the (leading) coefficients of 
Head(/) as / ranges over elements of B. By Noetherianity of A, we have J = (ci, . . . ,Cr)A for 
some Ci € A. Also, since r^diex is a well-partial-order, the final segment F = {Head(/) : / G 5} is 
finitely generated by T = {Head(6i), . . . ,Head(6|T|)} for some hj € B. Consider now the finite set, 

G := {ahj : 1 < i < r, 1 < J < |r[} C B. 

We claim that G is a subset of B fulfilling the requirements of the theorem statement. 

Let / G Br]A[mf. Then, Head(/ii) ^diex Head(/) for some hi £ {hi, . . . , with witness 
o"! € ©m (by Lemma 12. 9p . There are oi , . . . , A such that 

r 

/i := / - ^ OiCiaihi G B 

i=l 

has a strictly smaller (with respect to <diex) head than /. Continuing in this manner we can 
produce a sequence /i, /2, . . . of elements in B such that 

■ ■ ■ <diex Head(/2) <diex Head(/i) <diex Head(/). 

Since <diex is a well-ordering, it follows that fp = for some p G P which gives an expansion for / 
as in the statement of the theorem. □ 

Corollary 2.11. AlF]'' and A{F)^ are Noetherian A[&p]-modules. 

Remark 2.12. It turns out that Corollary 12.111 holds when A{F)^ is replaced by AS and A[(3p] 
by A[&] for any (3-set S with a nice ordering (this follows from the argument above). However, to 
prove Theorem 11.31 in the next section, we need the more refined statement found in Theorem l2.101 
which asks for witnesses a to (j2.ip having special properties. 
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3. Laurent chain stabilization 

In this short section, we prove that invariant chains of Laurent lattice ideals stabilize (this is 
Theorem 11.31 from the introduction) . The proof uses the order theory from the previous section 
and a few properties of lattice ideals. Some basic material on lattice and toric ideals can be found 
in \30\ Chapter 7] and [37J, respectively, and a more general reference for binomial ideals is |i6j. 

Let Q he a, finitely generated abelian group and let ai, . . . ,0^ be distinguished generators of Q. 
Let L denote the kernel of the surjective homomorphism Z"^ onto Q. The lattice ideal associated 
with L is the following ideal in ]K[zi, . . . , Zf^]: 

II = {z"" -z'' :u,v e n'^ with u-v£L). 

Here, we use the shorthand = z^^ ■ ■ ■ z^"^ for u = (ui, . . . , Ud) € Z*^. A toric ideal is the special 
case of a lattice ideal in which the group Q is torsion-free; in this case, the ideal II is also a prime 
ideal. 

Notice that if S* = {si, . . . , Sd} is a set with d elements, there is a natural isomorphism between 
Z"^ and the free Z-module ZS" with basis S given by: 

d 

(ai, . . . , ad) € Z"^ 1-^ OjSj € ZS. 

i=l 

Although simple, this identification will be useful for us below. 

Example 3.1. In the case S = (3)^ = {(1, 2), (1, 3), (2, 1), (2, 3), (3, 1), (3, 2)}, the integer vector 
(-2, 2, 1, 0, -1,0) G Z6 is also represented by -2 • (1, 2) + 2 • (1, 3) + (2, 1) - (3, 1) G Z(3)2. □ 

For simplicity of exposition, we focus our attention on lattice ideals in the polynomial rings 
TZji (equipped with the action of (3„) from (jl.3p . each of which has dn = indeterminates. Let 
Ln ^ Ln+i be an increasing sequence of subgroups of Z"'" C Z^"+i and let /„ := 7^,^ C (resp. 
In ^ ) be the corresponding lattice (resp. Laurent lattice) ideals. 

The basic idea in our proof of Theorem 11.31 is to view L = IJneP ^ Sp-invariant subgroup 
of the free abelian group Z[Ff = \J^^pZ[nf, which has free basis [Ff over Z. The set L has a 
finite generating set up to (3p-symmetry (using Theorem 12.101 and the fact that Z is Noetherian) , 
and these vectors are all contained in L^r for some integer A^. The remainder of the proof converts 
this fact back to the level of ideals. The complete details are as follows. 

Given an integer vector h E Z"^, we set G N*^ and /i_ € N'^ to be the nonegative and nonpositive 
part of h, respectively (so that h = — h-). The following is elementary. 

Lemma 3.2. Suppose that v,hi, . . . , hm G Z"' and set u = v + X^I^Li There exists a monomial 
z" £ K[zf\...,z^^] such that z^(z" -z^) G {z^^+ - z^^- : i = I, ■ ■ ■ ,rn)^^^±i^^^^^^±iy 

Proof. We shall induct on m, the base case being vacuously true. Consider the identity: 

(3.1) (z" -z^) = z^'-(z"-'^'" -z'')+z''~^"^-{z^"'+ -z^^-). 

As u' = u — hm has fewer terms, the proof follows by induction. □ 
Collecting these facts together, we can now prove the main result of this section. 
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Proof of Theorem ] l.iA The submodule L C Zp]'^ is finitely generated over Z[(3p] by Theorem l2.1UI 
Set H = G L) —G for a finite set of generators G L satisfying tlie property in Theorem 12.101 and 
let be such that H C Z"^^. Consider two vectors u,v e N"^™ such that n — f G Lm = Ln [m]'', 
with m > N. By assumption, the vector u — v is a Z-linear combination of ©^-permutations of 
elements in H. By Lemma 13.2^ it follows that z" — z" is a monomial multiple of an element in the 
ideal (of 7^^) generated by permutations (in 6m) of {z^+ — : h G H}. Thus, C {&mlN)fi± 
and the chain stabilizes. □ 

4. Stabilization of chains induced by monomials 

We now focus on the polynomial rings Rn from and the corresponding chains of toric ideals 
encountered in the statement of Theorem 11.41 

Definition 4.1. Let /c G P and / G . . . , y^]. For for each n > k, consider 

0n • Rn ^ Tji, («!,..., lift) ^ fi^ui :■■■■> tu^ )• 

Let In = ker (pn- The invariant chain Q Ik+i Q ■ ■ ■ is called the invariant chain of ideals induced 
by the polynomial f. 

The ideals in Definition 14. II appear in voting theory [TT], algebraic statistics [381 [251 [12l [17] , arid 
toric algebra [H [25] . When / is a monomial, each — ker <f>ji is a homogeneous toric ideal. The 
following was conjectured in 

Conjecture 4.2 ([!]). The chain of ideals induced by any monomial stabilizes. 

The authors of [4] verified the special case of Conjecture 14.21 when / is a square- free monomial. 
Underlying their proof is the fact that for every n > k, the ideals /„ are generated by quadratic 
binomials |37' Theorem 14.2]. Unfortunately, the corresponding statement is false when / is not 
square-free. Although a proof for the general conjecture is not known. Theorem 11.31 shows (albeit 
nonconstructively) that the Laurent versions of these chains stabilize. 

The main goal of this section is to provide an effective version of Theorem 11.31 for this situation 
that allows for explicit computation of generators (this is Theorem 11.41 from the introduction). In 
the next section, we describe this algorithm and give a reference to an implementation of it in 
software. We also explain another approach to finding these generators at the end of this section. 

Our running example throughout will be the case / = yfy2, and all computations were performed 
using Macaulay2 [18J and 4ti2 [Ij. If Iq is the chain of ideals induced by 2/^2/2; then Theorem 
11.31 guarantees stabilization of I^. Moreover, Theorem 11.41 provides a stabilization bound N = 
2 • deg(/) = 6. Using Algorithm [T] from Section [5l the following is a generating set for (below, 
we use a shorthand notation for indices; e.g., xie = g)): 

= {xiexl^X54,X65 - Xi4,Xi5x'^qX5Q, xj^X^^XisXQ^ - Xisxj^Xi^X'^Q, 

2 4 2 2 4 2 2 

■^16'^21"^45^65 ■^14'^15'^26' 2;i62;2]^ 2^342^65 2714X15X20X36, 

2 2 2 2 2 2 2 

XieX2iX3Q — X-^^^X2Q, XigX2iX32 — X12X13X26, 

XUX43-XMX34, a;i3X24 - 2;i4X23}. 

Therefore, the chain of Laurent ideals induced by yfy2 is generated by these 8 elements of G^ 
up to the action of the symmetric group. It is important to remark that these binomials are not 
generators of the original ideal Iq, nor of the chain /q. Moreover, this generating set is not smallest 
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possible, as shown in Section 14.21 where we study the combinatorial structure of this special case 
and find a generating set with only 4 elements for the Laurent chain I^. 

4.1. Proof of Theorem 11.41 First observe that the inclusion Rn gives us for every n> k 

an extension of (pn given by the homomorphism ipn ■ Rn ~^ '^n satisfying tpnixu) = (pnixu) and 
tpnix^^) = (pnixu)^^ for all u G (n)^. Notice that = kevipn and that we have the following 
commutative diagram: 

(4.1) Rn^ ^Rt 

Let a € be the exponent vector of a (non-constant) monomial / = = y"^ • • • y^*" , and 
consider 

An ■■= {cT(ai, . . . , afc, 0, . . . , 0)^ G : (J G 6„}. 

The set of column vectors An can be represented as an n x matrix with rows indexed by the 

indeterminates ti (for i = 1, . . . ,n) and columns indexed by the indeterminates Xw (for w G {n)^). 
The matrix An defines a semigroup homomorphism N" which lifts to the homomorphism 

The kernel In is generated by the set: 

jx'^ - : An{a) = An{h), a, h G n(^)^'} . 

For more details about toric ideals and their generating sets, see |37j . 

Example 4.3. Let A; = 2, n = 3, and a = (2, 1). The following represents the matrix ^3 associated 
to the homomorphism 03 defined by / = y'ly2- 

X12 xiz X21 X23 X31 X32 

ti 2 2 1 1 

t2 1 2 2 1 

t3 1 1 2 2 

The ideal I3 is generated by binomials: {a^isa^^i "^12^23; a;f3X2i — xf2a;3i, 3^212:31 — 3:12X32, X21X32 — 
xi2xl^, xi3a;23 - a;i2a;32, a;i3a;2ia;32 - a;i2X23a;3i, a;?3a;32 - xuxl^, X23XI1 - X13XI2, xl^X3i - a;2ix|2}- 

□ 

Next, we argue that it suffices to study those maps (pn ■ Rn — ^ Tn defined by an exponent 
vector a = {ai, . . . ,ak) G with gcd(a) := gcd(ai, . . . , at) = 1. To see this, suppose that 
gcd(a) = i > 1, and consider a' = ■ a. Let and (j)'^ be the homomorphisms given by 

<t>n{xw) = Ci • • 'Ct ^'^d (p'nixw) = tZ\ ■■■tZl, respectively. Note that (j)n{xw) = (</'n(^w))^ for all 
w G so if a, 6 G n(^)*^', then </.„(x«) = </>„(x^) ^ ^^(x'^)^ = c^'ni^^'Y ^ ^'^{^-) = </>Ux^) 

(as (f)'^{yi°') and <j}'^{y^) are monomials in T„); thus, x" — x^ G ker (pn if and only if x" — x^ G ker (p'^. 

Our first basic tool is a combinatorial lemma describing the Z-linear column span of An inside 
Z". For a G Z^, we set \a\ := Yli=i c^i- 
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Lemma 4.4. Let a = (ai,...,afc) £ with gcd(ai, . . . , a^) = 1. The integral span of the 
columns An {n > k) is: 

SpanzMn) = {/3 G : = mod \a\}. 

Proof. Let 2t = {/3 € : |/3[ = mod |a[}; clearly, Spanz(^„) C 21. By assumption, gcd(Q) = 1; 
thus, there are integers bi, . . . ,bk € Z with biai + • • • + 6^0^ = 1. 

For every j = 1, . . . ,n — 1 and i = 1, . . . , k, let aij be the transposition (i j) S S„, and consider 
the vector 

hj = bi{aija) H h bk{akja). 

Notice that hj is a vector whose jth. entry is 1 and hj G Span{A„}. Consider also the transposition 
Tj = {j n), and the vector 

f^'j = '^jhj = hiTjCTija) H h bk{Tjakja). 

For every i and j, the composition Tjaij is the transposition (in) € ©„; thus, h'j is obtained from 

/ij by changing the 1 from position j to position n. Naturally, h'- G Spana(^„). Let hj = hj — h'j G 

Span2(yl„). Notice that /ij is the vector with 1 in the jth position, —1 in the nth, and zeroes 
elsewhere. 

Now, let /3 = (/?!, . . . , /3„) G 21. By assumption, there exists g G Z such that |/3| = glaj. For every 
J = 1, . . . , n — 1, there is Vj G Z with /3j = qaj + rj. Set 

ra-l 

7 = ga + ^ Tj/ij G Spanz(^„). 

i=i 

It is easy to check that f3 = ^, and so /3 G Spanz(^„) as desired. □ 

Example 4.5. Consider a = (2, 1) and n = 3. Since gcd(2, 1) = 1, we can write 1 = (1)2 + (—1)1. 
The vectors /ii,/i2 £ Span^j^a} from the proof of Lemma HTil are precisely (1,0,-1)^, (0,1,-1)^. 
Therefore, the vector u = (1, 1, l)""" G 21 can be written as 

G Spanz(^3)- □ 

One immediate consequence of Lemma [4.41 is that the toric ideals in this section are not normal. 
This likely contributes to the difficulty of proving stabilization for chains induced by a non-square- 
free monomial. 

Corollary 4.6. Let be a non- square- free monomial in . . . , yj,]. For every n > \a\, the 

toric ideal In induced by the monomial y" is not normal. 

Recall from j37[ Proposition 13.5] that a toric ideal /_4 is normal if and only if pos{A) PI 
Spanz{A) = Spanfq{A), where pos{A) is the polyhedral cone defined by the columns of A. 

Proof. Let a G N'^ with gcd(a) = 1, and let r G &n be the cyclic permutation r = (12 - • • |a|). 
Realize a G Z" by a = (ai, . . . , a^, 0, . . . , 0)^ G Z". Consider the following identity: 

(1,...,1,0,...,0)^ = -^(a + ra + --- + Tl"l-2a + rl°l~ia). 

a 
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By construction z = (1, . . . , 1, 0, . . . ,0)^ € pos{An) and \z\ = \a\; thus, by Lemma HH] we see that 
z E pos{An) n Spanx{An)- However, since y° is not square-free, we must have z ^ SpanN(^„). □ 

Although not required for the proof of Theorem ll.41 the Smith normal form (SNF) of the matrices 
An can be easily computed from Lemma 14.41 For basic properties and algorithms involving the 
SNF over a principal ideal domain, we refer the reader to \22\ I39j. 

Corollary 4.7. Let a € such that gcd(a) = 1. For n > k consider the matrix 

An = {a{ai , . . . , afc, 0, . . . , 0)T G Z" : (J G ©„) . 
The Smith normal form for An is diag(l, . . . , 1, |a|). 

Proof. Use vectors hj from the proof of Lemma [4.4l to reduce the matrix An to SNF diag(l, . . . ,l,d), 
for some d G N. From Lemma 14.41 we know that Spa,nz{An) is 21 = {/3 G Z" : |/?| = mod \a\}. 
Since Z"/2l is a finitely generated Z-module, the fundamental decomposition theorem for modules 
[221 Theorem 7.8.2] implies that 

Z"/2l ^ Z/dZ. 

On the other hand, 21 is the kernel of the map Z" — )■ Z/|a|Z given by /3 i— > mod \a\; therefore, 

Z"/2t ^ Z/\a\Z, 

as Z-modules. Hence, Z/|a|Z = Z/dZ, which implies d = \a\. □ 
Let d = deg / = \a\ and set r = max{ai, . . . , Ofc}. Consider now 

(4.2) Bn := {(ai,... ,a„)^ G Z" : ai + • • • + a„ = d, < oi, . . . , < r}. 

There is a natural bijection between elements of Bn and multisubsets of [n] of cardinality d with 
at most r repetitions. Let r„ be the set of such multisubsets. Every a = (oi, . . . , a„) G is in 
bijection with a G Fn via: 

(4.3) a = {ai, . . . ,an) = {l''\2''\ . . . ,n''-}. 

Let Rn ■= K[Xr„]. When Bn is viewed as a matrix with rows indexed by ti (for i G [n]) 
and columns indexed by Xa (for a G F„), it defines a semigroup homomorphism that lifts to a 
homomorphism of K-algebras: 

4'n ■ Rn ^ ^71 • 

By definition, ^„ C Bn, and this inclusion gives an embedding r] : Rn ^ Rn- Also, extends the 
map (f)n in the sense that (j)n = 4>n° V- Therefore, we have the following commutative diagram: 



</>n 

Tn 

Example 4.8. Let n = 3 and a = (2,1). Then R^ = IK[2;i23, 2;ii2, 2:113, 2:122, a^23, 2:133, X233], and 
the following table represents the matrix B3 associated to the homomorphism (^3: 
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X123 


2;ii2 


a:ii3 


2^122 


2^223 


2;i33 


3^233 


h 


1 


2 


2 


1 





1 







1 


1 





2 


2 





1 




1 





1 





1 


2 


2 



□ 

We next note the following fact, easily derived using [37,, Theorem 14.2], as it provides a quadratic 
reduced Grobner basis for /„. These generators can be obtain from the quadratic generators of any 
Grobner basis for 

Lemma 4.9. The ideal C is generated by the quadratic binomials of any Grobner basis. 

In particular, since finite Grobner bases always exist, /„ has a finite set of quadratic binomials 
generating it. 

We now explain the key idea in our proof of Theorem 11.41 Since the map (pn extends (pn, we 
have In ^ In- Suppose that I„ = for some set Gn C and that we can find a K-algebra 
homomorphism fi making the following diagram commutative: 

(4.4) 




n 

Then, as is easily checked, fJ-{Gn) will be a generating set for I^. If, in addition, the Gn can 
themselves be finitely generated up to symmetry and fi is equivariant|3, then we have generated 
the whole Laurent chain I^ up to the symmetric group. As the proof of the following proposition 
explains, the existence of such a ;U is guaranteed by Lemma |4.4[ 

Proposition 4.10. Fix a = (ai, . . . , a^) € N" and let / = y" 7^ 1. For each n > k, there exists 
an equivariant K-algebra homomorphism /i : Rn — )• Rn that makes the diagram ()4.4p commute. 

Proof. Consider a multisubset a G F^. If € rj{Rn) ^ then define ^i{xa) '■= V'^i^a)- Assume 
Xd ^ rj{Rn). Since a is in bijection with a £ Bn as in (j4.3p . we have |o| = \a\. By Lemma 14.41 we 
can find integers B = {61, . . . , 6a/} C Z such that 

M 

a = y^^bjUj, 

i=l 



'The term equivariant for the map fi signifies that ii{ah) — aii{h) for any a € (Sn and h € R, 
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with M = {l)k\ and m G An- Let B+ = {k e B : bi e Z>o} and B_ = {k e B : b^ e Z<o}. 
Consider the fraction 



(4.5) c - ^"^^^^ 



Clearly q G R^, so we can define /u(xa) := q G i?^. 

Extend /i to i?„ by linearity. By construction, // makes the diagram ()4.4p commute since for 
a G r„, one can verify that ^'^(/"(a^a)) = </>n(2:a) = nr=i ■ ^ 

Remark 4.11. The above construction of /x is not necessarily unique as it depends on the repre- 
sentation of q. 

Example 4.12. Continuing from Example 14.51 we want to map X123 G i?3 to a fraction in that 
only involves the indeterminates of i?3 corresponding to those of R^,: 



123) 



XU3 



We also want this fraction to have the same image under as X123 has under ^3. Indeed, we have 
03(a;i23) = ht2t3 and 

, / . NX '03(a;il2a;33l) (/'3(2;i2)<?!>3(a;3l) (tft2)(iiil) , , , 
tp3[fJ'{xi23)) = —} ^ = -7 ^ = 77277 = tit2t3- 

'03(a;ii3) 03(a::i3) [tits) 

□ 

We are finally in position to prove Theorem II. 4[ 
Proof of Theorem \1.4\ Let /„ = ker(/)„; this ideal is generated by binomials of the form 

in which a U b ■ ■ ■ Uc = a' Li b' ■ ■ ■ Uc' as a union of multisets [S7\ Remark 14.1]. From Lemma |4.9| 



there is a finite generating set Qn of consisting of quadratic binomials. Let G„ be a finite set of 
generators for /„. Note that ri(In) Q In and so r]{Gn) C /„. For g G G„, we can write 

(4.6) T]{g) = ^ hpp, with hp £ Rn- 

We know G„ is a generating set for I^, but we give another generating set for in terms of Qfi . 
Applying the map ^ from Proposition 14. IOI to both sides of expression (j4.6p . we have 

P&Qn 

Moreover, //(p) G Ij^. It follows that = (^(p) : p G Qn) Since p G is a quadratic binomial, 

p = XaXi — Xa'Xy^ with o U 6 = o' U 6' as multisets. 

The cardinality of each of a, 5 is (i = |q|, and so the number of distinct numbers in aU 6 is at most 
2d. In particular, p.{p) G (©„/^)^± for n > 2d. Thus, stabilizes with bound N = 2d. □ 
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Example 4.13. Continuing with Example let g = X39X79 — X37X97 € Ig. Under the inclusion 
r] : Rq ^ Rq, we have r]{g) = X339a;779 — 2:337X799. From Lemma the ideal Ig is generated by 
quadratic binomials. We can write r]{g) in terms of those generators; in this case, 

riig) = {xssgXjjg - X379) - (X337X799 - xljg) £ Ig. 

Let pi := X337X779 — Xg^g and p2 ■= ^337X799 — X379. We have pi = aqi and p2 = crq2 for the 
following qi,q2 G h (actually /3 in this case) and o" = (1 3 9)(2 7) € ©9: 

qi = a;ii3a:223 - 2:123, q2 = 2:112X233 - 2;?23- 

Thus, fj-ipi) = cr/i(gi) and ^(^2) = 17^(92) since ^ is equivariant and fj.{qi) and fJ-{q2) generate g up 
to symmetry: 

(™2 ™2 \ / ^2 ™2 
X12X22. 2"^ ~ ^ ( ^12^32 2^ 
2^13 / V ^13 

□ 

4.2. Toric ideals induced by y\y2- Theorem 11.41 provides evidence that chains of ideals induced 
by monomials stabilize. The simplest (unknown) case is when / = yiy2- Here, we present an 
explicit computation of the generators for the corresponding Laurent chain that is different from 
Algorithm [H We hope to illustrate some of the complexity of the general problem and also to 
elaborate on other approaches for tackling Conjecture 14.21 

For n > 2, let /„ be the toric ideal induced by the monomial yiy2- Let An G be the 

matrix that defines the semigroup homomorphism (j)^ such that In = ker (recall Definition 14. ip . 
For example, when n = 5 we have 

/22221111000000000000\ 
10002000222111000000 
01000200100200221100 
00100020010020102021 
VOOOIOOO2OOIOO2OIO212/ 



^5 



When the columns of An are ordered lexicographically, a basis for \eii[An) as a Z-module can 
be described as follows: 



(4.7) kerz(A 



where X^-n is the (d — n) x (d — n) identity matrix and /C^ is a matrix whose structure we now 
describe. Let cv G Z"'"^ be the row vector whose entries are all equal to r. Then, 

= ( -^1 £^2 -^3 ^a) 1 

in which 
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For instance, when n = 5, the integer kernel of A5 has the following Z-basis 



2 


2 


2 


2 


2 


2 


3 


3 


3 


4 


4 


4 


4 


4 


] 


-2 
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For each n, the elements of kerz(>t„) are Z-linear combinations of the columns of the matrix 
(j4.7p . For each i = 1, ... ,4, we can realize the columns of Ci as the first column of Ci after 
applying a permutation a G 6„ to it. For instance, when n = 5, the first column of Ci is the vector 
(2, —2, 0, 0, —1, 1,0,..., 0)^ G Ij"^^, which corresponds to the binomial xfga^si — Xi^X2i. If we apply 
the transposition (3 5) G 65 to this element, we get x^2^5i ~ ^15^21, whose corresponding integer 
vector is precisely the third column of Ci; namely, (2, 0, 0, —2, —1, 0, 0, 1, 0, ... , 0)^ G Z^°. 

In general, for every n and for i = 1,2,3, the transposition (3j) with A < j < n applied to the 
binomial corresponding to the first column of Ci will be equal to the binomial whose support corre- 
sponds to the (j— 2)-th column of Ci. For £4, instead of transpositions, we use those permutations 
that send the pair (3,4) to («, j) for 3 < i 7^ j < n to write those binomials corresponding to the 
columns of £4 in terms of the first column of C^^. For instance, the binomial x\2X^4, — 3:14x^3X21 has 
support the first column of £4. When we apply the permutation (3 4 5) G ©5 to this binomial, we 
get X12X45 — X15X14X21, which has support the 5th column of £4. Consider the set 

TT± r 2 2 2 23 224 22i 

H = 1x12X31 - X13X21, X12X23 - 2:13X21, X12X32 - X13X21, X12X34 - 3:13X14X211 

of binomials corresponding to the first column of each Ci. The action of 65 on produces 
generators for the Laurent ideal corresponding to the toric ideal /s, by Lemma |3.2[ In general 
for n > 5, the action of ©n on produces generators for the Laurent ideal I^. We thus obtain a 
generating set for the chain that depends only on the description of keTz{An) and is independent 
from the methods used in the proof of Theorem 11.41 

Unfortunately, we could not generalize this technique to other cases as the combinatorics that 
describe kerz(^n) in general becomes more complicated. We also remark that the set fails to 
be a generating set for the (non-Laurent) chain of ideals induced by yfy2- 
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5. Algorithms 

The proof of Theorem 11.41 suggests an algorithm to find the generators of a chain of Laurent 
toric ideals induced by a monomial y". We stated the existence of such an algorithm in Theorem 
11.71 from the introduction. In this section we describe this algorithm and argue its correctness. A 
full implementation in Macaulay2 [18] can be found in |24] . 



Algorithm 1 [Theorem 1 1.7| 
Input: Exponent vector a G N'^ 

Output: Generators for the chain of Laurent ideals defined by y° up to symmetry 
1: d := 2\a\ 

2: Compute the matrix Bd (j4.2p 

3: Compute the Grobner basis G of the toric ideal Ig^ 
4: for all g £ Q do 

5: for all indeterminates in g do 
6: if not indexed by a permutation of a then 

7: g = replace Xw in g by the monomial quotient ^{x^) 

8: end if 
9: end for 
10: end for 

11: Remove redundant generators from Q 
12: return Q 



Given an exponent vector a E N'^, the algorithm computes a set of generators for the chain 
of Laurent ideals defined by up to the action of the symmetric group. In the first steps, it 
considers all the integer partitions oi d = 2|a| with parts at most maxct := max{ai, . . . ,an}) and 
then constructs the matrix Bd by taking all the permutations of such partitions. 

In step m the algorithm constructs the toric ideal Id that corresponds to the matrix Bd and 
computes its Grobner basis Q (with respect to any term order). This Grobner basis computation 
is the most expensive step for large ideals. We decided to use the Macaulay2 package FourTiTwo, 
which invokes one of the fastest routines, 4ti2, specializing in computing Grobner bases for toric 
ideals [T]. 

Step [TT] removes the redundant generators from Q. Using Lemma 14. 9( we start by removing 
all the non-quadratic generators from Q. We then remove the symmetric orbit of each of the 
remaining generators. To illustrate how drastically the number of generators is decreased after 
this step, consider once more the running example of Section [H When y" = y\y2-, the Laurent 
toric chain has a stabilization bound at n = 6; for this value of n, the toric ideal /g C Rq has 
270 minimal generators. When we lift to the ideal /g C Rq, we obtain 849 minimal generators, 
but only 13 modulo the action of the symmetric group. From those, we find that 11 generate the 
corresponding Laurent ideal modulo the symmetric group. But after clearing denominators and 
common monomial factors, we found that only 8 from those 11 (exactly those 8 that are presented 
in the introduction of Section H]) form a generating set of the Laurent ideal modulo the action 
of the symmetric group. Since the number of generators increase when passing to the ring one 
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way to improve speed on this orbit removal step is to remove the orbit after Step [3] and again in 
Step El 

The core of our algorithm is Step U\ where we turn Lemma 14.41 and Proposition 14.101 into a 
computational tool. We unfold this step in Algorithm [2] below. This algorithm expresses every 
element of the column span of Bd as a linear combination of the column span of the matrix Ad 
using the construction found in the proof of Lemma 14. 4^ and detailed in Algorithm [3] below. This 
integer decomposition is then used to create the map fj, in (14. ID . 



Algorithm 2 Construction of the map /i 

Input: indeterminate x^j, exponent vector a G N*^ with gcd(Q) = 1 
Output: monomial quotient fi{xw) 

1: Set V := {aa : a G ©„} 

2: if w then 

3: Write w = bivi + • • • + bj-Vr with 6j G Z and Vi G V for all i G [r] (Lemma 14. 4p 

4; index := {i G [r] : hi > 0} and index := {i £ [r] : bi < 0} 

5: numerator := 1 and denominator := 1 

6; for all i G indexB^ do 

7: numerator = numerator ■ x!** 

8: end for 

9; for all i G indexB^ do 
10; denominator = denominator ■ x~^^ 
11: end for 

12: return numerator /denominator 
13: else 

14: return x^ 
15: end if 



We remark that in the union computation in Step 7 of Algorithm [Sj we add coefficients of 
matching pairs as in {oo-, Va} U {b^, v^} = {aa + ba,Va}- 

6. Open problems and conjectures 

Stabilization of chains of ideals is unexpected and important for applications. However, the 
problem of deciding whether a chain is stable under the action of a group seems difficult, even for 
the special case of the symmetric group. In this section, we present some conjectures based on 
computational evidence. We focus first when the ideals comprising the chain are toric ideals as 
they tend to have rich combinatorial structure; later, we turn to a more general setting and close 
with some problems that develop this topic further. 

Motivated by the study of bounds on the Castelnuovo-Mumford regularity in algebraic geometry, 
Bayer and Mumford introduced in [5j the degree- complexity of a homogeneous ideal I with respect 
to a term order ^ as the maximal degree in a reduced Grobner basis of I, and this is the largest 
degree of a minimal generator of in-<(/). In our context, degree-complexity is important because it 
is closely related to stabilization of chains of ideals. For instance, in the proof of Theorem 11.41 we 
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Algorithm 3 Integer decomposition of /3 in terms of a 

Input: Integer vector /3 G Z*^, exponent vector a G N*^ with gcd(a) = 1 and |a| dividing |/3| 
Output: List {{a(j,Vcr} : cr € ©„} such that /3 = X](tg6„ ^o-'^cr, where fg- = o" • a, and Oo- € Z 
1: Write 1 = h\a\ + • • • + fefcCKfc, where a = (ai, . . . , a^) and bi £ Z 

2: g:= |/3|/|al 

3: L := {{q,Va}}, where = (qi, . . . , ctfe, 0, . . . , 0) G Z" 
4; for j from 1 to n — 1 do 
5: if /3j — Qfj • g 7^ then 
6: for i from 1 to /c do 

7: L = LU {{{f3j - aj ■ q)bi,Va^^}, {-bi{/3j - aj ■ q),Vrja,j}}, where aij and tj are as in the 

proof of Lemma 14.41 

8: end for 

9: end if 
10: end for 
11: return L 



exploited the fact that the ideal /„ is binomial and has degree-complexity 2 for every n. On the 
other hand, if the ideals comprising a chain induced by a monomial do not have a degree-complexity 
bound, then stabilization is unlikely. 

We pose the following problems based on our observations in Table [6] (computed using our 
software pi]). 

Conjecture 6.1. Let a = (ai,a2) with gcd(ai,Q;2) = 1 (suppose qi > 02)- The degree- complexity 
of In is of the form 2a\ — 02 for all (but possible finitely many) ideals In in the chain of (non- 
Laurent) toric ideals induced by the monomial y". 

Problem 6.2. Let a G with gcd(a) = 1; is the degree- complexity of In constant as n —t- 00? 

Recall the set i3„ from (j4.2p . Using the fact that Iq^ is generated by quadratics we show in this 
paper that for An ^ Bn, the chain of ideals 7^,^ has a corresponding Laurent chain that is 
stable under the action of 6p. On the other hand. Conjecture 14.21 makes the stronger claim that the 
chain stabilizes. While it is difficult to find subsets C„ C Bn for which the chain of ideals Ic„ 
is stable under the action of Sp, one might get some indications by solving the following problem. 

Problem 6.3. Find combinatorially defined subsets Cn ^ i3„ such that the toric ideal Ic„ has 
constant degree- complexity as n grows. 

This problem is of particular interest in algebraic statistics. For instance, in (21[ Conjecture 7.3], 
it is conjectured that for any T > 3 and a fixed 5 > 3, the toric ideals of the homogeneous Markov 
chain model on S states are generated by polynomials of degree at most S — 1. This is an instance 
of Problem 16.31 as for each T > 3, the design matrix of such a model is precisely a subset of the 
matrix Bn for n = T. 

The early stabilization of structured chains appears to be common. It would be interesting to 
construct examples of chains with nontrivial lower bounds on stabilization. 
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Table 1. Degree-comp 



exity of the toric ideal In defined by y*^ 



Problem 6.4. Let f{d) be an increasing function / : P — )• P. Find a family of invariant chains 

{(d) 1 ^ 
lo f (over TZf or Rp) which have stabilization bound at least f{d). 

More specifically, we ask whether a linear lower bound holds for the chains in Theorem 11.41 and 
their Laurent counterparts. 

Jo > from Theorem \1.4\ must 

have stabilization bounds at least f{d) = Cd. 
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